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Abstract. 

We study the topological susceptibility, X; in QCD with two quark flavours using the 
lattice field configurations recently produced by the UKQCD collaboration. We find 
clear evidence for the expected suppression at small quark mass, and examine the 
variation of x with this mass. The resulting estimate of the pion decay constant, 
= 105 ± 5 IlIq MeV, is consistent with the experimental value of ~ 93 MeV. We 
compare x to the large- iV^ prediction and find consistency over a large range of quark 
masses. The benefits of the non-perturbative action improvement scheme and the 
matched lattice spacings between simulation ensembles are discussed. We compare our 
results with other studies and suggest a reason why such a quark mass dependence has 
not been previously seen. 
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1 Introduction 



The ability to access the non-perturbative sectors, and to vary parameters fixed in 
Nature has made lattice Monte Carlo simulation a valuable tool for investigating QCD 
and related theories, and especially so as far as the role of topological excitations is 
concerned. 

In gluodynamics (the pure gauge or "quenched" theory) lattice calculations of the 
continuum topological susceptibility now appear to be relatively free of the systematic 
errors arising from the discretisation, the finite volumes and the various measurement 
algorithms employed. Attempts to measure the microscopic topological structure of 
the vacuum are also well advanced (for a recent review, see IQ). 

The inclusion of sea quarks in ("dynamical") lattice simulations, even at the rela- 
tively large quark masses currently employed, is numerically extremely expensive, and 
can only be done for lattices with relatively few sites (typically 16^.32). To avoid 
significant finite volume contamination of the results, the lattice must be relatively 
coarse, with a spacing a ~ 0.1 fm. Compared to quenched lattice studies at least, 
this is a significant fraction of the mean instanton radius, and has so far precluded a 
robust, detailed study of the local topological features of the vacuum in the presence of 
sea quarks. The topological susceptibility, on the other hand, may be calculated with 
some confidence and provides one of the first opportunities to test some of the more 
interesting predictions for QCD. Indeed, it is in these measurements that we find some 
of the most striking evidence for the presence of the sea quarks (or, alternatively, for a 
strong quenching effect) in the lattice simulations. 

In this talk we study the topological susceptibility on ensembles of SU(3) lattice 
gauge fields that have been recently generated by the UKQCD collaboration using a 
QCD lattice action with Nf = 2 degenerate fiavours of (fully dynamical) sea quarks. 
These ensembles have been produced with two notable features. The first is the use of 
an improved action, such that leading order lattice discretisation effects are expected 
to depend quadratically, rather than linearly, in the lattice spacing (just as in gluody- 
namics). In addition, the action parameters have been chosen to maintain a relatively 
constant lattice spacing, particularly for the larger values of the quark mass. 

We find these features allow us to see the first clear evidence for the expected sup- 
pression of the topological susceptibility in the chiral limit, despite our relatively large 
quark masses. From this behaviour we can directly estimate the pion decay constant 
without needing to know the lattice operator renormalisation factors that arise in more 
conventional calculations. 

The structure of this talk is as follows. In Section |^ we discuss the measurement 
of the topological susceptibility and its expected behaviour near the chiral limit, and 
in the limit of a large number of colours, Nc. In Section ^ we describe the UKQCD 
ensembles, and the lattice measurements of the topological susceptibility over a range 
of sea quark masses. We fit these with various ansatze motivated by the previous 
section. We compare our findings with other recent studies in Section ^. Finally, we 
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provide a summary in Section |[ 

An early version of the results presented here can be found in p[, While the results 
presented in this talk are still preliminary, we do not expect any significant changes to 
our conclusions. The final complete version, including estimates of the r( mass, will 
appear in 0. 



2 The topological susceptibility 

In four-dimensional Euclidean space-time, SU(3) field configurations can be separated 
into classes , and moving between different classes is not possible by a smooth deforma- 
tion of the fields. The classes are characterised by an integer-valued winding number. 
This Pontryagin index, or topological charge, Q, can be measured by integrating the 
local topological charge density 

Q{x) = (1) 

over all space time 

Q = ^ [d'x-Q{x). (2) 

The topological susceptibility is the expectation value of the squared charge, normalised 
by the volume 

X-f^- (3) 

Sea quarks in the vacuum induce an instanton-anti-instanton attraction which becomes 
stronger as the quark masses, m^, . . . , decrease towards zero (the "chiral" limit) 
and this leads to a suppression of the topological charge and susceptibility to leading 
order in the quark mass 0] 

X = S(— + — (4) 

where 

S = - lim lim {0\i>'^\0) (5) 

is the chiral condensate (see for a recent discussion). Here we have assumed 
{0\ijjilj\0) = (OImmIO) = (0|(i(i|0) and we neglect the contribution of heavier quarks. 
PCAC theory relates this to the pion decay constant and mass as 

f^ml = (m„ + md)S + 0(mj) (6) 

and we may combine these for Nj degenerate light flavours to obtain 

X = ^ + Oimt) (7) 
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in a convention where the experimental value of /tt ~ 93 MeV. This relation should 
hold in the limit flm^V ^ 1. We anticipate our results here to say that even on our 
most chiral lattices the LHS is of order 20 and so this bound is well satisfied. Thus a 
calculation of x as a function of m^r will allow us to obtain a value of f^^. This method 
has an advantage over more conventional calculations in that it does not require us to 
know difficult-to-calculate lattice operator renormalisation constants. 

As iriq and m^r increase away from zero we expect higher order terms to check the 
rate of increase of the topological susceptibility so that, as m^, m,r ^ C)0, it approaches 
the quenched value, x'^"- fact, as we shall see below, the values of x that we obtain 
are not very much smaller than x^". So there is the danger of a substantial systematic 
error in simply applying Eqn. ^ at our smallest values of uIt^ in order to estimate /^r. 
To estimate this error it would be useful to have some understanding of how x behaves 
over the whole range of m,. This is the question to which we now turn. 

There are two quite different reasons why x might not be much smaller than x^i". 
The obvious first possibility is that irtq is 'large'. The second possibility is more subtle: 
ruq may be 'small' but QCD may be close to its large- A^^c limit 0. Because fermion 
effects are non-leading in powers of Nc-, we expect x ~^ X^" fo^' any fixed value of m^, 
however small, as the number of colours Nc oo. There are phenomenological reasons 
P, for believing that QCD is 'close' to Nc = oo, and so this is not idle speculation. 
Moreover in the case of D=2+l gauge theories it has been shown Q that even SU(2) 
is close to SU(cx)). Cruder calculations in four dimensions indicate that the same is 
true there. For example the SU(2) and SU(3) quenched susceptibilities are very similar 
[m . Since our lighter quark masses straddle the strange quark mass it is not obvious if 
we should regard them as being large or small. We shall therefore take seriously both 
of the possibilities described above. 

We start by assuming the quark mass is small but that we are close to the large- A'^c 
limit. In this limit, the topological susceptibility is known ||^ to vary as 



where x°°j /oo are the quantities at leading order in A^,,- In the chiral limit, at fixed N^ 
this reproduces Eqn. in the large- A'c limit, at fixed m^, it plateaus to the quenched 
susceptibility, x°°y because oc A'c. The corrections to Eqn. || are of higher order in 
and/or lower order in Nc- 
We now consider the alternative possibility: that rUq is not small, that higher order 
corrections to x will be important for most of the values of rUq at which we perform 
calculations, and that we therefore need an expression for x thai interpolates between 
rUq = and rUq = oo. Clearly we cannot hope to derive such an expression, so we will 
simply choose one that we can plausibly argue is approximately correct. The form we 
choose is 
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3 



The coefficients have been chosen so that this reproduces Eqn. |^ when m^r —>■ and 
X X^" + ^{^/''^t) when 771^^ oo. Thus this interpolation formula possesses the 
correct limits and it approaches those limits with power-like corrections, as we would 
expect. (The precise form of the correction as m.,^ —>■ oo will also depend on the physical 
quantity that is used to set the overall mass scale; we shall not attempt to determine 
that here.) 

We shall use the expressions in Eqns. | and || to analyse the dependence of 
our calculated values of x and to obtain a value of /^r together with an estimate of the 
systematic error on that value. In addition the comparison with Eqn. ^ can provide us 
with some evidence for whether QCD is close to its large- A'^c limit or not. 



3 Lattice measurements 

We have calculated x fo^i' complete ensembles of dynamical configurations, pro- 
duced by the UKQCD collaboration |10|, as well as on one which is still in progress. 



The SU(3) gauge fields are governed by the Wilson plaquette action, with "clover" 
improved Wilson fermions. The improvement is non-perturbative, with Csw chosen 
to render the leading order discretisation errors quadratic (rather than linear) in the 
lattice spacing, a. 

The theory has two coupling constants. In pure gluodynamics the gauge coupling, 
/3, controls the lattice spacing, with larger values reducing a as we move towards the 
critical value at /3 = oo. In simulations with dynamical fermions it has the same role 
for a fixed fermion coupling, k. The latter controls the quark mass, with k ^ Kc 
from below corresponding to the massless limit. In dynamical simulations, however, 
the fermion coupling also affects the lattice spacing, which will become larger as k is 
reduced (and hence increased) at fixed /?. 

The three least chiral UKQCD ensembles (by which we mean largest m^r/mp) are 
calculated using couplings (/?, k) = (5.29,0.1340), (5.26,0.1345) and (5.20,0.1350) (in 
order of decreasing sea quark mass). By appropriately decreasing /3 as k is increased, 
the couplings are "matched" to maintain a constant lattice spacing |Tl], [1^] (which is 
'equivalent' to /3 ~ 5.9 in gluodynamics with a Wilson action) whilst approaching the 
chiral limit. Discretisation and finite volume effects should thus be of similar sizes on 
these lattices. The fourth ensemble at (5.20,0.1355) and the preliminary results that 
we present for (5.20,0.13565) have lower quark masses, but are at a reduced lattice 
spacing. To have maintained the matched lattice spacing here would have required a 
reduction in /?. As the non-perturbative value Csw(/5) is not known below P = 5.20 the 
matching has had to be sacrificed for the continued removal of the 0{a) effects. As the 
lattices are all L > 1.5 fm, we believe that the minor reduction in our lattice volume 
should not lead to significant finite volume contamination of results. 

Four-dimensional lattice theories are scale free, and the dimensionless lattice quan- 
tities must be cast in physical units through the use of a known scale. For this work. 
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we use the Sommer scale [|12| both to define the lattice spacing for the matching pro- 
cedure, and to set the scale. The measured value of fo on each ensemble corresponds 
to the same physical value of vq = 0.49 fm. (rg is the dimensionless lattice value of tq 
in lattice units i.e. fg = ro(a)/a. We use the same notation for other quantities.) As 
we are in the scaling window of the theory, we can then use the naive dimensions of 
the various operators to relate lattice and physical quantities, e.g. f^x = f^oX + C'(a^), 
where we have incorporated the expected non-perturbative removal of the corrections 
linear in the lattice spacing. 

Further details of the parameters and scale determination are (to be) given in [^. 
Measurements were made on ensembles of 400-800 configurations of size L^T = 16^.32, 
separated by ten hybrid Monte Carlo trajectories. Correlations in the data were man- 
aged through jack-knife binning of the data, using bin sizes large enough that neigh- 
bouring bin averages may be regarded as uncorrelated. 

We begin, however, with a discussion of lattice operators and results in the quenched 
theory. 

3.1 Lattice operators and x^" 

The simplest lattice topological charge density operator is 

Qi^) = ^£fiuaT^rU^^{n)U^r{x) (10) 

where U^y{n) denotes the product of SU(3) link variables around a given plaquette. 
We use a refiection-symmetrised version and form 

Q = (11) 

X - 1^ (12) 

with L^T the lattice volume. 

In general, Q will not give an integer-valued topological charge due to finite lattice 
spacing effects. There are at least three sources of these. First is the breaking of scale 
invariance by the lattice which leads to the smallest instantons having a suppressed 
action (at least with the Wilson action) and a topological charge less than unity (at least 
with the operator in Eqn. p!0| ). We do not address this problem in this study, although 
attempts can be made to correct for it |13[, but simply accept this as part of the overall 



0{a ) error. In addition to this, the underlying topological signal on the lattice is 
distorted by the presence of large amounts of UV noise on the scale of the lattice spacing 



14| , and by a multiplicative renormalisation factor |T5[ that is unity in the continuum, 
but otherwise suppresses the observed charge. Various solutions to these problems 
exist [ffl]. In this study we opt for the 'cooling' approach. Cooling explicitly erases 
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the ultraviolet fluctuations so that the perturbative lattice renormalisation factors for 
the topological charge and susceptibility are driven to their trivial continuum values, 
leaving O(a^) corrections that may be absorbed into all the other lattice corrections 
of this order. We cool by moving through the lattice in a "staggered" fashion, cooling 
each link by minimising the Wilson gauge action applied to each of the three SU(2) 
subgroups in the link element in turn. (The Wilson gauge action is the most local, and 
thus the most efficient at removing short distance fluctuations whilst preserving the 
long range correlations in the flelds.) Carrying out this procedure once on every link 
constitutes a cooling sweep (or "cool"). The violation of the instanton scale invariance 
on the lattice, with a Wilson action, is such that an isolated instanton cooled in this 
way will slowly shrink, and will eventually disappear when its core size is of the order 
of a lattice spacing, leading to a corresponding jump in the topological charge. Such 
events can, of course, be detected by monitoring Q as a function of the number of 
cooling sweeps, ric- Instanton-anti-instanton pairs may also annihilate, but this has no 
net effect on Q. All this does, however, motivate us to perform the minimum number 
of cools necessary to obtain an estimate of Q that is stable with further increasing Uc 
(subject to the above). In our calibration studies we found {Q){nc) and (Q^)(nc) to be 
stable within statistical errors for 5, and the results presented here are for = 10. 

We also remark that on a lattice one obviously loses instantons with sizes p < 0{a). 
Since the (pure gauge) instanton density decreases as this would appear to induce a 
negligible 0{a^) error in the susceptibility. It can be numerically substantial, however, 
for the coarse lattices often used in dynamical simulations. 

In general, then, we expect the topological charge and susceptibility to be suppressed 
at non-zero lattice spacing. In gluodynamics with the Wilson action this can be fltted 
well by just a leading order, and negative, correction term 

f^X = 0.072 - 0.208/f2 (13) 

in the window of lattice spacings that is comparable to those used in current dynamical 
simulations. We shall refer to this formula in what follows. 

3.2 Sea quark effects in the topological susceptibility 

In Table |l] we give our estimates of the topological susceptibility. We now convert 
these raw lattice values into physical units using tq as the physical scale. In Fig. |I| 
we plot the resulting values of versus a similarly scaled pseudoscalar meson mass 
(calculated, of course, with valence quarks that are degenerate with those in the sea, 
i.e. Kvaience = '^sea)- For each poiut we also plot the corresponding value of the quenched 
topological susceptibility, calculated at the same lattice spacing using the interpolation 
formula Eqn. |T^. It is useful to note at this point the small variation in the quenched 
values over the range of data points simulated. This is a consequence of the UKQCD 
strategy of matching fg. Without this, if the calculations were performed at flxed (3, the 
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lattice spacing would become increasingly coarse with increasing m^, and the reduction 
in the quenched susceptibility would be much more pronounced over this range of k. 
We shall return to this important point when we discuss other work in a later section. 

Comparing the dynamical and quenched values, the effects of the sea quarks are 
clear. Whilst the (5.29,0.1340) point is consistent with the quenched value, moving to 
smaller rriq (oc m\) the topological susceptibility is increasingly suppressed. 

We can make this observation more quantitative by attempting to fit our values of 
with the expected functional form in Eqn. |^. But we must first be clear whether 
this fit is justified, and what exactly we are extrapolating in. Eqn. ^ is strictly a chiral 
expansion that describes the behaviour for small sea quark masses in the continuum 
and relates it to /^r. We expect that it is also applicable for a set of lattice data points 
evaluated at the same lattice spacing, except that now the decay constant, will 
be the one appropriate to that lattice spacing. Whilst most of our data points are 
evaluated on a trajectory of constant lattice spacing in the parameter space, the two 
most chiral values are not. The lattice spacings there are slightly finer. If fo/7r varied 
significantly with a over this range of a, it would not be clear how to perform a consis- 
tent chiral extrapolation through the data points. The non-perturbative improvement 
of the action, however, removes the leading order lattice spacing dependence and the 
residual corrections in this range of lattice spacings appear to be small, at least in 



measurements of hadron spectroscopy [O, |T^. Assuming the variation in the decay 



constants is similarly small over what, it should be recalled, is not a large variation 
in a, we may proceed to attempt a common chiral extrapolation to the data. For 
this purpose it is useful to redisplay the data in Fig. |^, where the leading order chiral 
behaviour would then be a horizontal line, 

^ = CO (14) 

and including the leading correction gives a general straight line 

^ = Co + ci(fom^)^. (15) 

In each case the intercept is related to the decay constant by cq = (fo/7r)^/4. We now 
follow a standard fitting procedure, first using the most chiral points, then systemati- 
cally adding the less chiral points until the fit becomes unacceptably bad. The larger 
the number of points one can add in this way, the more evidence one has for the fitted 
form and the more confident one is that the systematic errors, associated with the 
neglected higher order corrections, are small. The results of performing such fits are 
shown in Table ^ and those using the two and four most chiral points respectively are 



plotted on Fig. ^. We see from the Table that the fits using Eqn. |15| show much greater 



stability and these are the ones that will provide our eventual best estimate for /^r. 
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We should comment briefly on the determination of the fitting parameter errors. In 
performing all but the fiat fit we must contend with the data having (small) errors on 
the abscissa in addition to the ordinate. In order to estimate their affect on the fitting 
parameters, we first perform fits to the data assuming that the abscissa data take their 
central values. Identical fits are then made using the central values plus one, and then 
minus one standard deviation. The spread of the fit parameters obtained provides what 
is probably a crude over-estimate of this error (given there is some correlation between 
the ordinal and abscissal uncertainties) but is sufficient to show that it is minor. We 
show this spread as a second error, and for estimates of the decay constant we add it 
in quadrature to the other fit parameter error. 

It is remarkable that we can obtain stable fits to most of our data using just the 
first correction term in Eqn. |15|. Nonetheless, as we can see in Fig. |l], our values of 
the susceptibility are not very much smaller than the = oo quenched value and 
we need to have some estimate of the possible systematic errors that may arise from 
neglecting the higher order corrections that will eventually check the rise in x. As 
discussed earlier we shall do so by exploring two possibilities. One is that the reason 
why X is close to x^^ is not that is 'large' but rather that A'^c = 3 is large. Then the 
values of x should follow the form in Eqn. A second possibility is simply that our 
values of nig are indeed large. In that case we have argued that the functional form 
Eqn. ^ should be a reasonable representation of the true mass dependence. We now 
perform both types of fit in turn. 

We begin with the first possibility, and therefore fit the data with the following 
ansatz 

ro X= : p-—^^ 16 

where we expect C3 = fo^x"^^ up to 0{1/N^) corrections. To test this we fit up to 
six data points. The first five are measured in the dynamical simulations. The final 
quantity is the quenched susceptibility at a lattice spacing equivalent to those of the 
three least chiral, "matched" ensembles. The last is assigned an arbitary large value 
(fo'^Tr)^ = 1000 whose exact value has no effect on the fits. 
We also expect, from the Maclaurin chiral expansion of Eqn. |T6|, 

ro'^X = Co ■ {rorrinY - cl/cs ■ (rom^)^ + cj^/cg ■ {rQin^f + 0{{rQm^f) (17) 

that Co is related to the decay constant as before, Cq = (fo/7r)^/4. We present the 
results of the fits in Table ^ 

We turn now to fits based on the functional form in Eqn. |^. We therefore use the 
ansatz 

fo'x = ^(rom.)^ tan- (^T^^) (18) 
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where once again we expect C3 = ro x'^" ^ind from the expansion 

rVx = coirom^f - (^)' ■ (^) ■ (romj^ + O{{fom^f) (19) 

we expect Cq = (ro/7r)^/4. Note that in contrast to Eqn. Eqn. ^ has no term that 
is cubic in rUg and the rise will remain approximately quadratic for a greater range in 
(fom^)^. That this need be no bad thing is suggested by the relatively large range over 
which we could fit Eqn. [1^. And indeed we see from the fits listed in Table ^ that this 



form fits our data quite well. 

Typical examples of the fits from Eqn. 16 and Eqn. 18 are shown in Figs. and ^ 



The similarity of the two functions is apparent. In Table ^ we use the fit parameters to 
construct the first three expansion coefficients in the Maclaurin series for the various 
fit functions, describing the chiral behaviour of x- The fits are consistent with one 
another. 

The fitted asymptote of the susceptibility at large rriq is given by C3. We see from 
Table ^ that these are broadly consistent with the quenched value, and our large 
statistical errors do not currently allow us to resolve any 0{1/N^) deviation from this. 

As an aside we ask what happens if we cast aside some of our theoretical expectations 
and ask how strong is the evidence from our data that (a) the dependence is on 
rather than on some other power, and (b) the susceptibility really does go to zero as 
—>■ 0? To answer the first question we perform fits of the kind in Eqn. |1^ but 
replacing (fom^)^ by (ro^Tr)^- We find, using all six values of x, that c = 1.64(22) (8); 
a value consistent with c = 2. The x^/d-o.f. is poorer, however, than for the fit with 
a power fixed to 2 (possible as there is one fewer d.o.f.) suggesting that the data does 
not warrant the use of such an extra parameter. As for the second question, we add a 
constant c to Eqn. |18| and find c = 0.0096(78)(20). Again this is consistent with our 
theoretical expectation; and again the x^/d.o.f. is worse. (See Table ^ for details of 
the above two fits.) 

Given the consistency of our description of the small regime from our measure- 
ments it is reasonable to use the values of Cq to estimate the pion decay constant, f.^. 
This is done in units of tq in Tables and ^. We use the chiral fit of Eqn. |15| over the 
largest acceptable range to provide us with our best estimate and its statistical error. 
We then use the fits with other functional forms to provide us with the systematic 
error. This produces an estimate 

foA = 0.261 ±0.013 1°;°^^ (20) 

where the first error is statistical and the second is systematic. This is of course no 
more than our best estimate of the value of /^r corresponding to our lattice spacing of 
a ~ 0.1 fm. This value will contain corresponding lattice spacing corrections and these 
must be estimated before making a serious comparison with the experimental value. 
Our use of the non-perturbative Cs„ should have eliminated the leading 0{a) errors. 
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however, and (a/ro)^ is sufficiently small that it is plausible to expect these lattice 
corrections to be smaller than our other errors. If we use the phenomenological value 



Tq = 0.49 fm then we obtain from Eqn. ^ the value 

U = 105 ±5 t\o MeV (21) 

which is reasonably close to the experimental value ~ 93 MeV. In a comparison 
with other physical quantities is expected to provide us with a quantitative estimate 
of the magnitude of the O(a^) errors. 



4 Comparison with other studies 

During the course of this work, there have appeared some other studies of the topo- 
logical susceptibility in lattice QCD; in particular one by the CP-PACS collaboration 
|T8| , and one by the Pisa group |T^. Both calculations find no significant decrease of 



the susceptibility with decreasing quark (or pion) mass when everything is expressed 
in physical, rather than lattice, units. This appears to contradict our findings and is 
clearly something we need to address. 

Our suggestion why these other studies have seen no decrease in x is as follows. They 
differ from our study in having been performed at fixed f3. That implies that the lattice 
spacing a decreases as ruq is decreased. In typical current calculations this variation in 
a is substantial. (See for example Fig. 4 of pO|.) At the smallest values of nig the lattice 
spacing cannot be allowed to be too fine, because the total spatial volume must remain 
adequately large. This implies that at the larger values of rriq the lattice spacing is 
quite coarse. Over such a range of lattice spacings the topological susceptibility in the 
pure gauge theory typically shows a large variation. Since for coarser a more instantons 
(those with p < 0{a)) are excluded, and more of those remaining are narrow in lattice 
units (with a correspondingly suppressed lattice topological charge) we expect that this 
variation is quite general, and not a special feature of the pure gauge theory with a 
Wilson action. In lattice QCD we therefore expect two simultaneous effects in x as we 
decrease at fixed /3. First, because of the C(a^) lattice corrections just discussed, x 
will, like x^", tend to increase. Second, it will tend to decrease because of the physical 
quark mass dependence. In the range of quark masses covered in current calculations 
this latter decrease is not very large (as we have seen in our work) and we suggest that 
the two effects may largely compensate each other so as to produce a susceptibility 
that shows very little variation with m^; in contrast to the ratio x/x'^^ which does. 

To illustrate this consider the fixed-/? calculation in [^]. The range of quark masses 
covered in that work corresponds to (rom^r)^ decreasing from about 6.5 to about 3.0. 
Simultaneously a/ro decreases from about 0.437 to about 0.274. Over this range of 
1/fo = a/ro the pure gauge susceptibility increases by almost a factor of two, as we see 
using Eqn. |l^. Clearly this is large enough to compensate for the hoped-for near-chiral 
variation of the susceptibility. 
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To make this a little more quantitative we consider the different calculations in turn, 



beginning with that in [|I9|. The Pisa group uses the Wilson gauge action, and four 
flavours of staggered fermions. Simulations are carried out at a fixed value of the gauge 
coupling for a range of bare quark masses between rhg = 0.01 and 0.05. Expressed in 
units of the string tension, a, this corresponds to a change in the quark mass by 
a factor of approximately three. Over the corresponding range of lattice spacings the 
susceptibility in the pure gauge theory with Wilson action varies by a factor of 1.44(19), 
when determined using the same "field theoretic + smearing" method (for a review, see 
0). Such an increase would be large enough to largely compensate for the hoped-for 
near chiral decrease of the susceptibility (see for example the variation of f^x between 
(fom,r)^ — 6.0 and ~ 2.0 in Fig. 

The CP-PACS study is carried out using an improved (Iwasaki) gauge action with 
two flavours of tadpole-improved Sheikoslami-Wohlert clover fermions. Combining 
the data from [|l^, one finds over the range of quark masses considered a variation 
in the square of the pion mass (in units of the string tension) of approximately 3, 
similar to that of the Pisa study. Over this range there is also little variation in the 
topological susceptibility when expressed in units of the string tension. The string 
tensions considered would correspond, in a pure gauge Wilson action study, to a range 
in P of 5.64 — 5.71, and in that range the topological susceptibihty, x^^, varies by 
a factor of around 1.5. If we assume that there is a similar effect with the Iwasaki 
action, then over this range the dynamical x appears to suffer a suppression of about 
a factor of 1.5 relative to the pure gauge theory. As before, not all the data is near the 
chiral limit, and over a corresponding range in the UKQCD data we see a suppression 
of around a factor of 2. CP-PACS in fact possess datasets at four values of /3 and 
they have calculated the topological susceptibility on the dataset corresponding to the 
second smallest value. At their higher values of (3 the effect we have been describing 
here will be much weaker and we believe that calculations performed on those two 
datasets will reveal a quark mass dependence comparable to the one we find. 

The arguments presented in this section are by no means definitive. We do not, 
for instance, know the discretisation effects on the topological susceptibility in lattice 
gluodynamics with the Iwasaki action. Neither do we know know how the (different) 
fermionic actions alter these discretisation effects. What the heuristic analysis given 
here aims to show is that these results are not necessarily in contradiction to the 
theoretical expectations and to the results presented in this talk. 



5 Conclusions 

We have calculated the topological susceptibility in lattice QCD with two light quark 
flavours, using lattice field configurations, recently generated by UKQCD, in which the 
lattice spacing is approximately constant as the quark mass is varied. We find that 
there is clear evidence of the effects of sea quarks in suppressing x- 
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We discuss this behaviour in the context of chiral and large Nc expansions, and find 
good agreement with the functional forms expected there. We are not able to make a 
stronger statement about how close QCD is to its large limit, owing to the relatively 
large statistical errors on our calculated values, particularly at larger quark masses. 

The consistent leading order chiral behaviour from our various fitting ansatze allows 
us to make an estimate for the pion decay constant, /j^ = 105 ± 5 MeV, for the 
lattice spacing of a ~ 0.1 fm. (Here the first error is statistical and the second has to 
do with the chiral extrapolation.) Since we use a lattice fermion action in which the 
leading 0{a) discretisation errors have been removed, and (quenched) hadron masses 
show little residual lattice spacing dependence, we might expect that this value is close 
to its continuum limit. In any case we note that it is in agreement with the experimental 
value, ~ 93 MeV. 

We have commented on some recent studies of the topological susceptibility in lattice 
QCD, which do not see the kind of mass dependence that we claim to have observed. 
We have suggested that this might arise from the fact that in these fixed-/? calculations 
the value of a decreases with decreasing and the consequent variation in those 
C(a^) lattice corrections that the susceptibility (plausibly) shares with the pure gauge 
theory might be large enough to compensate for the expected near-chiral variation. 
Further analysis of this and the corresponding quenched theories in this light would be 
welcome. 
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X X 10^ 


(5.20, 0.13565) 
(5.20, 0.13550) 
(5.20, 0.13500) 
(5.26, 0.13450) 
(5.29, 0.13400) 


3.18 (64) 
4.84 (57) 
7.90 (42) 
8.86 (54) 
12.9 (1.9) 



Table 1: The topological susceptibility in lattice units. 



Fit 


iVfit 


Co 


Cl 


xVd.o.f. 




Eqn. 14 


2 


0.0140 (16) 




0.805 


0.237 (14) 


Eqn. |14 


3 


0.0112 (6) 




2.202 


0.212 (6) 


Eqn. 14 


4 


0.0091 (4) 




9.008 




Eqn. |15] 


3 


0.0176 (35) (4) 


-0.0018 (10) (1) 


0.964 


0.265 (27) 


Eqn. PI 


4 


0.0170 (16) (1) 


-0.0016 (4) (0) 


0.502 


0.261 (13) 


Eqn. 15 


5 


0.0147 (14) (1) 


-0.0011 (3) (0) 


2.965 


0.242 (12) 


Table 2: Fits to the A^gt most chiral points 


of{rlx)/rhl. 


Fit 


iVfit 


Co 


C3 


xVd.o.f. 


fofn 


Eqn. 


3 


0.0208 (87) (12) 


0.0844 (427) (35) 


1.013 


0.288 (61) 


Eqn. 16 


4 


0.0272 (85) (18) 


0.0632 (114) (6) 


0.895 


0.329 (53) 


Eqn. |1^ 


5 


0.0233 (66) (10) 


0.0717 (147) (3) 


1.847 


0.305 (44) 


Eqn. H 


6 


0.0279 (32) (11) 


0.0631 (5) (0) 


1.489 


0.334 (21) 


Eqn. 


3 


0.0186 (53) (7) 


0.0576 (175) (6) 


0.990 


0.273 (40) 


Eqn. H 


4 


0.0209 (42) (7) 


0.0506 (55) (5) 


0.682 


0.289 (30) 


Eqn. |T| 


5 


0.0189 (36) (5) 


0.0550 (69) (6) 


1.929 


0.275 (27) 


Eqn. H 


6 


0.0164 (14) (7) 


0.0629 (5) (0) 


1.567 


0.256 (12) 



Table 3: Fits to the Nat most chiral points of (^qx). 



Fit 


iVfit 


const. 


O((rom.)2) 




Eqn. 15| 


4 


0.0170 (17) 


-0.0016 (4) 




Eqn. |1| 
Eqn. |1| 


5 
6 


0.0233 (67) 
0.0279 (34) 


-0.0031 (19) 
-0.0123 (30) 


0.0025 (24) 
0.0055 (19) 


Eqn. 18] 
Eqn. |18 


5 
6 


0.0189 (37) 
0.0164 (16) 


-0.0026 (11) 
-0.0017 (4) 







Table 4: Chiral expansion terms of fitted functions. 
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iVfit 


Co 


C3 


c 


c 


xVd.o.f. 


6 


0.0112 (40) (5) 


0.0533 (78) (17) 


fixed to 2 


0.0096 (78) (20) 


1.825 


6 


0.0206 (26) (16) 


0.0631 (80) (19) 


1.61(22) (8) 


fixed to 


1.798 



Table 5: Fits to the Nfn most chiral values of (^qx); for the power c ofm^, and for 
the chiral intercept, c. These are as described in the text. 
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0.08 



0.06 



0.04 



0.02 



0.00 




o dynamical x 
A interpolated x"*"' 
- fit to Eqn. 16 
fit to Eqn. 18 



Figure 1: The measured topological susceptibility, with interpolated quenched points at 
the same Tq and fits independent of the quenched points. 
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0.03 



0.02 



0.01 




■ fit to Eqn. 14 
A fit to Eqn. 15 
Vfit to Eqn. 16 
O fit to Eqn. 18 
* pliysical value 



0.00 ' ' ' ' ' ' ' ' ' ' ' 

2 4 6 8 

Figure 2: Fits to the measured topological susceptibility, independent of the quenched 
points. 
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